Production of intense beams of mass-selected water cluster ions and 
theoretical study of atom-water interactions 
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(Dated: March 30, 2009) 

The influences of water molecules surrounding biological molecules during irradiation with heavy 
particles (atoms, ions) are currently a major subject in radiation science on a molecular level. In 
order to elucidate the underlying complex reaction mechanisms we have initiated a joint experimental 
and theoretical investigation with the aim to make direct comparisons between experimental and 
theoretical results. As a first step, studies of collisions of a water molecule with a neutral projectile 
(C atom) at high velocities (> O.la.u.), and with a charged projectile (proton) at low velocities (< 
O.la.u.) have been studied within the microscopic framework. In particular, time-dependent density 
functional theory (TDDFT) was applied to the valence electrons and coupled non-adiabatically to 
Molecular dynamics (MD) for ionic cores. Complementary experimental developments have been 
carried out to study projectile interactions with accelerated (< lOkeV) and mass-selected cluster 
ions. The first size distributions of protonated water cluster ions H+(H20)n(n=2-39) produced using 
this new apparatus are presented. 

PACS numbers: xxx 



I. INTRODUCTION 



While the use of ionizing radiation is well- 
established, notably in therapies and analytical 
techniques, each new development opens a field of 
investigation around the possible dangers to our 
health and the cnvironmcntlv? . Concurrently, be- 
sides applications and risk evaluations, the effect 
of ionizing radiation in biomolccular nanosystcms 
is emerging as a major area of research, both on 
a fundamental level and as a source for experi- 
mental and technical innovations. The irradiation 



*Dept. Physics and Astronomy, The Open University 
(OU), Walton hall, Milton Keynes, MK76AA, UK 
t Universite Hassan II-Mohammedia, Facult des Sciences 
Ben M'Sik (LPMC), B.P.7955 Ben M'Sik, Casablanca, Mo- 



of biomolccular nanosystems in the gas phase un- 
der single collision conditions represents a signif- 
icant new subject in radiation science. From an 
experimental standpoint, it has great potential for 
the development of new analytical techniques and 
synthesis techniques. Complementary theoretical 
progress can provide new descriptions of radia- 
tive energy transfer mechanisms on the molecular 
scale, opening new perspectives for the elucidation 
of the radiation dose in living systems. 

Our principle goal is to investigate microscopic 
mechanisms responsible for irradiation effects in 
biological molecules with a particular focus on 
the influence of water molecules in the environ- 
ment of these molecules. This is achieved here in 
a joined experimental and theoretical effort with 
the aim to make direct comparisons between ex- 
perimental and theoretical results. The originality 
of this work lies in the ability to quantify the role 
of a biomolcculc's immediate environment in pro- 
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ton radiation-induced processes through the pre- 
cise control of the number of associated water 
molecules. This is a key aspect from the funda- 
mental physics point of view and it is also cru- 
cial to allow realistic and detailed comparisons 
between experiments and theory. The path to- 
wards such an ambitious goal is long and due to 
the many elementary processes involved one needs 
to validate the various experimental and theoreti- 
cal aspects. First of all, we consider here the sim- 
plest case of a single water molecule interacting 
with a projectile, which constitutes a prerequisite 
for any further developments. The next step in- 
volving an assemble of water molecules around a 
biomolecule raises no major theoretical difficulty 
because the number of molecules can be efficiently 
controlled. Monitoring this number experimen- 
tally constitutes a key issue for further experimen- 
tal developments and for comparisons with theory. 

In previous experiments at the Institut de 
Physique Nucleaire de Lyon, the ionization of gas- 
phase water molecules by protons^ and neutral 
hydrogen^ impact has been studied in the veloc- 
ity range of 0.9-2.5 a. u. (units of the Bohr veloc- 
ity), coinciding with the Bragg peak maximum for 
energy deposition by an ion beam in an absorb- 
ing medium. Mass-analyzed 1120"*" and fragment 
ions of the water molecules were detected in co- 
incidence with the projectile following ionization 
in single collision conditions. The determination 
of the charge state of the projectile after the colli- 
sion enabled ionization processes to be separated 
on the basis of charge transfer between the target 
molecule and the projectile. Thus it was possi- 
ble to measure branching ratios and absolute cross 
sections for water ionization with and without the 
capture of an electron by the incident proton. 

The present experimental challenge is to study 
interactions between water molecules in a clus- 
ter following proton impact induced excita- 
tion/ionization (20-150 keV)of water prior to dis- 
sociation. Accordingly a new experimental system 
has been developed to investigate the proton irra- 
diation of mass-selected protonated water clusters. 
Many different techniques have been used success- 
fully to generate ensembles of water cluster ions in- 
cluding electrospray ionization^i^, electron impact 
ionization^, corona discharge ionization^, chemical 
ionization** and electrospray droplet impact^iiiS. 
The novelty of the present experimental system 
lies in the acceleration (up to 10 kV) and mass- 
selection of molecular and cluster ions of water 



prior to 0.9-2.5 a.u. collisions with protons. 

The microscopic description of radiation induced 
processes requires an explicit dynamical account 
of electronic degrees of freedom which respond 
first in the present collisions. Moreover, it is nec- 
essary to treat electrons in a non-adiabatic way 
and to allow for ionization and/or electron trans- 
port. This invalidates most calculations based on 
the Born-Oppenheimer (BO) approximation ex- 
cept in some specific situations. Indeed, depend- 
ing on the characteristics of the ionizing projec- 
tile (charge, velocity), one can treat the problem 
in a simplified manner by decoupling electronic 
and ionic dynamics. A typical example is the 
case of high velocity charged projectiles in which 
ions can be safely considered as frozen and the dy- 
namics reduced to the electronic response, at least 
for short times. Another example is the case of 
low velocity neutral projectiles for which a ground 
state (Born-Oppenheimer, BO) treatment is ac- 
ceptable. To the best of our knowledge, low ve- 
locity charged and high velocity neutral projec- 
tiles can nevertheless not be treated by the usu- 
ally available calculationsii, and certainly not in 
the framework of a unique theoretical approacbi^. 
Fortunately, in the case of cluster dynamics, sev- 
eral approachesi^im^ii^iii have been developed 
which consider this question of coupled electronic 
and ionic dynamics in relation to irradiation by 
intense laser fields. A coupling of the optically 
active spot to a large environment can be added 
in a hierarchical approac h^^i^^ . In the present 
work we have adapted the non adiabatic approach 
of Calvayrac et al.— to the case of an organic 
molecule and applied it to realistic irradiation sce- 
narios. The method of Calvayrac et ali^ contains 
as limiting cases pure electron dynamics and BO 
dynamics (Car-Parinello dynamics) and can thus 
describe collisions with high-velocity charged or 
low- velocity neutral projectiles. Moreover, at vari- 
ance with currently available approaches, it en- 
ables complementing cases to be described such 
as low-velocity charged and high-velocity neutral 
projectiles. Indeed such a non adiabatic approach 
places no restriction on the velocity or charge state 
of the projectile and thus offers an unified picture 
of many possible irradiation scenarios. 
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chamber and detected through a 1 mm coUima- 
tor by a channehron electron muhipher operated 
in a counting mode. Further experimental devel- 
opments are currently being carried out in order 
to intersect the cluster ion beam with an intense 
beam of 0.9-2.5 a.u. protons or a jet of atomic or 
molecular gas (He, N2...). 



FIG. 1: Schematic view of the experiment 



II. EXPERIMENTAL SET-UP 



The present experimental system is shown 
schematically in Figure [TJ All measurements are 
taken with a two-sector-field mass spectrometer of 
EB geometry combined with an electron-impact 
ion source. Neutral clusters are produced by ex- 
panding water vapor from a stagnation chamber 
at a temperature of around lOOC and a pressure 
of approximately 1 bar through a pin-hole nozzle 
(20 fj.m diameter) into vacuum. In this isentropic 
expansion the temperature drops rapidly with in- 
creasing distance from the nozzle leading to super- 
saturation of the water vapor and subsequent clus- 
tering. The clusters are ionized by electrons of 
50 eV (experimentally available range from about 
0-50 eV) and the pressure in this region is main- 
tained to 2x10"'* mbar by a 12001/s oil diffusion 
pump. The resulting positive ions and cluster ions 
(charge state q) are immediately extracted from 
the center of the ion source by a weak penetrating 
electric field. The extracted ions pass through a 
skimmer and are accelerated to (6xq) kV. After 
the ionization chamber the ions pass a skimmer. 
An Einzel lens and four pairs of deflection plates 
are used to direct the beam to the focal point of a 
hemispherical electric sector field with maximum 
transmission. The electric sector field serves to 
select ions according to their kinetic energy, thus 
increasing the energy resolution of the beam. The 
KE-resolved ion beam then passes through a mag- 
netic sector field for ion selection according to mo- 
mentum. Finally, the mass-selected water clus- 
ter ions H+(H20)n are focused into the collision 



III. THEORETICAL METHODS 



In order to perform microscopic simulations of dy- 
namical processes, we employ time-dependent den- 
sity functional theory (TDDFT) for the electrons 
combined with classical molecular dynamics (MD) 
for the ionic cores. The TDDFT is used at the level 
of the local-density approximation (LDA) together 
with an average-density self-interaction correction 
(ADSIC) to achieve the correct ionization poten- 
tials (IP), a feature which is crucial to describe 
electron emission correctly. Wavefunctions and 
fields are represented on a 3D coordinate-space 
grid of dimensions 72*72*64. Absorbing boundary 
conditions are used to remove outgoing electrons. 
Thus the total number of electrons N = N{t) de- 
creases in time. The number of escaped electrons 
A^csc = N{t = 0) — N{t) is a measure of average ion- 
ization. The details of this method are presented 
elsewhere (e.g.^i^o and for ADSIG^i). 

The coupling between electrons and ionic cores is 
achieved by Goedecker-typo22, non-local pseudo- 
potentials (PsP). The original parameterizations 
employ different Gaussian widths for each material 
and contribution. This hampers numerical preci- 
sion in coordinate space grids. We have refitted 
new PsP parameters for the smallest width the 
same value of 0.412 ap in all terms (local and non- 
local) and for all elements involved in the study. 
Further details will be given in a forthcoming pub- 
lication. 
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FIG. 2: Spectrum of mass-selected water cluster ions 
at 6 keV energy obtained by scanning the magnetic 
sector field (see figure [TJ. 



size. It is worth noting that the intensity of the 
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FIG. 3: Spectrum of mass-selected water cluster ions 
at 6 keV energy in a log scale 



IV. FIRST RESULTS 



A. Experimental production of size resolved 
protonated water clusters 



Figure 2 shows a typical mass spectrum obtained 
for the water cluster beam by scanning the mag- 
netic sector field. This spectrum was recorded for 
an acceleration voltage of 6 kV and for a water va- 
por pressure and temperature of 1.1 bar and 100 
C in the stagnation chamber, respectively. The 
water monomer ion peak is not shown in the spec- 
trum as its intensity is too high leading to satu- 
ration of the channcltron detector. The observed 
series of peaks spaced by 18 amu / unit charge is 
attributed to protonated water clusters compris- 
ing n molecules (n = 2-39). Water cluster ion in- 
tensities decrease exponentially (red dashed line 
in Fig. 2) as the number of water molecules in- 
creases. The other peaks correspond to compo- 
nents arising from the residual gas in the cluster 
source. The insert in figure 2 shows the higher- 
mass peaks (n=18-39) in greater detail. These 
peaks correspond to H+(H20)n (n = 18-39) cluster 
ions. The present size distribution of cluster ions 
is rather similar to that observed by Castleman 
and co-workers2.. The present results are plotted 
in Figure 3 also on a logarithmic scale. In this case 
a clear drop in intensity can be seen for ion pro- 
duction above cluster size n=21. This is consistent 
with previous results^ supporting a magic number 
(a particularly stable configuration) at this cluster 



mass-analysed cluster ion beams achieved with the 
present set-up is rather high. For instance, the in- 
tensity of the protonated water trimer beam mea- 
sured through the 1 mm-coUimated channeltron is 
5000 counts/s. In a next step these mass selected 
cluster ion beams will be crossed by an intense pro- 
ton beam and the resulting collision events will be 
analyzed using a molecular imaging detection sys- 
tem. 



B. Theoretical studies of the irradiation of 
water molecules 



As discussed in section U the ultimate aim of our 
combined theoretical and experimental project is 
a systematic investigation of collision induced pro- 
cesses in H+(H2 0)n cluster cations. To achieve 
this ambitious goal, it is necessary to validate the 
employed methods at intermediate steps. As the 
key building block of H+(H20)n clusters, the neu- 
tral H2O molecule provides the starting point for 
our theoretical investigation. Systematic studies 
of small H+(Il20)n (n = 1-3) will be presented 
elsewhere. The first models developed here will de- 
scribe the collision of an atom or ion with an H2O 
molecule. As obscrvablcs, we study the ionization 
and the motion of the various atoms in the sys- 
tem. In particular we have explored the two limit- 
ing cases of low- velocity neutral and high- velocity 
charged projectiles for which one can reduce the 
dynamics to Born-Oppenheimer or pure electronic 
motion, respectively. We have checked that our 
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non-adiabatic approach enables us to recover these 
hmitinff expected. In the present stud- 

ies, we find that non-adiabatic effects do indeed 
remain negligible which a posteriori validates pre- 
vious investigations. We thus focus on the cases 
of low-velocity charged and high-velocity neutral 
projectiles for which an accurate coupling between 
ionic and electronic dynamics becomes compul- 
sory. We illustrate the capabilities of our approach 
on such cases. 



1. Neutral atom impact upon H2O at high-velocity 



Fig.|3]shows results for a collision of a fast (v=0.1- 
0.2 a.u.) neutral C atom with an H2O molecule. 

The water molecule lies in the scattering plane. 
The upper left panel shows the time evolution of 
ionization for impact velocity w=:0.1 a.u.=20 ao/fs 
and varying impact parameter b. The sudden 
jump to A'csc ~ 4 is caused by the C atom 
leaving the numerical box and carrying its four 
active electrons away. The net ionization of the 
water molecule is the difference A'esc-4. There is 
negligible ionization (and excitation) of the water 
molecule for the largest h (3.15 ao). For h = 2.21 
ao, however, we see a loss of about 0.2 charge 
units whereas for the smallest b (1.87 ap) one full 
electron is lost. The reason for this becomes clear 
when we consider the right upper panel showing 
the trajectories in the scattering plane for the 
smallest impact parameter. The C atom travels 
from right to left through the box with negligible 
deviations from a straight line, first triggering 
ionization by about 0.2 charge units as shown 
in the left-hand panel. At the point of closest 
impact, the C atom hits one of the H atoms and 
transfers a sufficiently large momentum that the 
H atom dissociates from the molecule and leaves 
the box taking the attached 0.8 electrons away. 

The lower right panel shows the trajectories 
for the impact parameter. There is again a strong 
impact on the H atom closest to the projectile. 
As for the smaller impact parameter, there is 
still a significant ionization of about 0.2 charge 
units. However, the transferred momentum is 
sufficiently low that the atom remains bound in 
the molecule. The H atom rotates around the 
O atom and pushes the other H atom into the 
same rotation such that the structure of molecule 



is essentially conserved. Thus atom impact has 
triggered here just a strong rotational motion. 
Finally, the lower left panel shows ionization 
for medium impact parameter h = 2.21 ag and 
two different velocities (0.1 a.u. and 0.2 a.u.). 
Ionization increases with velocity although the 
change is surprisingly small. We have here two 
counteracting influences. The faster atom exposes 
the molecule higher frequency components so it 
can enhance ionization. On the other hand, the 
time of interaction becomes much shorter leading 
to a reduced reaction yield. The combination 
of these effects leads to a rather small velocity 
dependence in this range. 



2. Proton impact upon H2 O at low velocity 



The calculations treating charged projectile 
collisions with H2O at a small impact parameter 
(6 = 0.81 ao) and at various low velocities are 
shown in Fig. [5l The II2O molecule is placed with 
the O atom in the scattering plane and the two 
H atoms facing out of the plane. The interaction 
between the charged projectile and the molecule 
has a larger range than the neutral projectile case 
due to ion charge coupling with the large dipole 
moment of the water molecule. As a consequence, 
the projectile trajectories experience bending 
depending on the projectile velocity. The upper 
right panel shows the lowest velocity case and 
it is clear that the bending angle is very high. 
The water molecule as a whole is accelerated 
to a translational motion combined with some 
rotation of the two H atoms relative to the O 
atom (compare the initial and final configurations 
shown in the inserts). The lower right panel 
represents a collision at 0.05 a.u.. The bending 
of the projectile trajectory and the effect on the 
motion of the II2O molecule are both markedly 
smaller than observed for v=0.025 a.u.. Again, 
the molecule is accelerated to translational motion 
with some small rotation, but this time in the 
other direction. The fastest collision is shown in 
the lower left panel. The much shorter interaction 
time reduces bending and impact on the molecule. 
The left upper panel shows ionization for the 
various impact velocities. The slowest collision, 
although having the largest effect on molecular 
motion, produces the least electron emission. 
Ionization increases with increasing velocity and 
seems to level off near the largest velocity. Thus, 



6 



HgO + C 




v=0.1 a.u. 

— b=3.15a, 

— b=2.21 a, 

— b=1.87a. 



5 10 
time [fs] 



15 



20 




rr 



b=2.21 ag 

— v=0.1 a.u. 

— v=0.2a.u. 



5 10 15 
time [fs] 



20 



15 
10 
5 

""o 

™, 
>, 

-5 
-10 
-15 



b=2.21 ao 


final 


v=0.1a.u. 




^^^^ 






• C 




• H 




• 



-15 -10 



-5 
x[ao] 



10 15 



FIG. 4: Modeled collision between a water molecule and a C atom at various high velocities and impact parameters 
h. Left panels : number of escaped electrons as a function of time. Right panels: ionic trajectories in the x — y 
plane; the arrows indicate the time evolution; the insert in the top right panel shows the initial configuration of 
H2O for all calculations, while that in the bottom right panel presents the final configuration only for the case 
V = 0.1 a.u. and b = 2.21 ao. 



the enhancement by higher frequency components 
dominates over reduced interaction time in that 
(low) velocity range. 

The theoretical test cases have demonstrated 
the feasibility of detailed dynamical simulations of 
vifater systems excited by various projectiles. We 
have seen a significant difference in the dynamics 
triggered by neutral (C atom) versus charged 
(H+ ion) projectiles. There is sizeable impact 
ionization. The effects on the molecular motion 
as a whole are large and depend sensitively on the 
impact conditions. 



V. CONCLUSION AND OUTLOOK 



In this paper we have demonstrated the abil- 
ity to produce intense beams of accelerated and 
size-selected water cluster ions in our newly con- 
structed cluster/atom collision apparatus. The ac- 
companying theoretical studies have demonstrated 
our ability to explore the dynamics of water sys- 
tems irradiated by various projectiles. The next 
steps of the experimental investigation will focus 
on the effects of irradiation of size selected wa- 
ter cluster ions. The parallel theoretical investiga- 
tions are currently being extended to consider ir- 
radiation of small singly charged water clusters. In 
particular, comparisons will be drawn between our 
experimental and theoretic studies of size-selected 
water cluster ion collisions with high- velocity pro- 
tons. Attachment of water molecules to a molecule 
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FIG. 5: Modeled collision between a water molecule and a H"*" ion at low velocities (v=0. 015-0.1 a.u.) for a 
fixed impact parameter b — 0.81 ao. Top left panel: number of escaped electrons as a function of time. Other 
panels: ionic trajectories in the x-z plane. The arrows indicate the time evolution; only in the top right panel, 
the arrows on the water molecule are drawn at same instant. The insert in the bottom left panel shows the initial 
configuration of H2O for all the calculations, while the other inserts present the final configuration for v = 0.05 
a.u. (bottom right) and v — 0.025 a.u. (top right). 



of biological interest will constitute a further step 
both from the experimental and theoretical sides. 
The experimental capability to control the number 
of attached water molecules is essential for a direct 
comparison with theoretical calculations. This op- 
portunity will be exploited in the studies to come. 
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